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Abstract 

It is shown that the contribution from double parton scattering to the inclusive double heavy meson yield is quite 
comparable with the usually considered mechanism of their production at the LHC energy. For some pairs of heavy 
flavored quarks in the final state the double parton scattering will be a dominant mode of their production. 
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In the last years it has become obvious that multiple 
parton interactions play an important role in hadron- 
hadron collisions at high energies and are one of the 
most common, yet poorly understood [1], phenomenon 
at the LHC. The presence of such multiple parton in- 
teractions in high-energy hadronic collisions has been 
convinciiigly demonstrated by the AFS [2], UA2 Isfc, 
CDF P, Ist], and DO |6] Collaborations, using events 
with the four-jets and y + 3-jets final states, thus pro- 
viding new and complementary information on the pro- 
ton structure. The possibility of observing two separate 
hard collisions was proposed long ago. Early theoreti- 
cal investigations were carried out in the framework of 
the parton model I?! S 01 with subsequent extension to 



literature (see, for instance, 

13, H [il s im in in 




iijisllil 



references therein). 

A greater rate of events containing multiple hard in- 
teractions is anticipated at the LHC with respect to the 
experiments mentioned above due to the much higher 
luminosity and greater energy of the LHC. Moreover 
the products from multiple interactions will represent 
an important background [12, 13, 14J to signals from 
the Higgs and other interesting processes and certain 
types of multiple interactions will have distinctive sig- 
nature 1 15, 16, 17, 18] facilitating a detailed investiga- 
tion of these processes experimentally. 



The main purpose of this letter is to bring attention 
to another important processes: the production of heavy 
meson pairs through double parton scattering which is 
definitely not taken into consideration in the current the- 
oretical estimations (3^ 31 ]. Here, however, one should 
mention quite a recent paper 1I32I1 . in which the con- 
tribution from the double parton scattering to y/i/r-pair 
production has been discussed for the first time for the 
condition of the LHCb experiment. 

Let us recall that, with only the assumption of fac- 
torization of the two hard parton processes A and B, 
the inclusive cross section of a double parton scattering 
process in a hadron collision is written in the following 
form 



"dps 



i,j,k,l 



I 



r,7(xi,x2;bi,b2;ei 



■ Ql) 



xdiixi , x'l , e^o-f/fe, x'2, Ql) 
xru(x\ , x'2, bi - b, b2 - b; Q^, Ql) 
xdx\dx2dx^dx2cf'bicf'b2cf'b, (1) 

where b is the usual impact parameter, that is, the 
distance between the centers of colliding (e.g., the 
beam and the target) hadrons in transverse plane. 
r,j(xi, X2; bi, b2; Q\, QI) are the double parton distribu- 
tion functions, depending on the longitudinal momen- 
tum fractions xi and X2 and on the transverse positions 
bi and bz of the two partons undergoing the hard pro- 
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cesses A and B at the scales Qi and Q2; tr^, and are 
the parton-level subprocess cross sections. The factor 
m/2 is a consequence of the symmetry with respect to 
the interacting parton species / and j: m - 1 if A = B, 
and m = 2 otherwise. 

It is typically taken that the double parton distribution 
functions may be decomposed in terms of the longitudi- 
nal and transverse components as follows: 



rij(xuX2-M,h2;QiQl) 
^D',^(xuX2;QlQl)f(hi)f(h2), 



(2) 



where /(bi) is supposed to be an universal function for 
all kind of partons with its normalization fixed as 



f(hi)f(bi-b)d^bid^b ^ I T(h)d^b^l, (3) 



and r(b) - J/(bi)/(bi - h)d^bi is the overlap func- 
tion. 

If one makes a further assumption that the longitudi- 
nal component D'i^(xi,X2', Q\, Q^) reduces to a product 
of two independent one parton distributions, 

£>;/(xi ,x2;Ql Ql) = Diixi ; Q])dI(x2; Q\\ (4) 

the cross section of the double parton scattering can be 
expressed in a simple form 



DPS 



(t" rr" 
2 o-„(r 



(5) 
(6) 



In this representation and at the factorization of longi- 
tudinal and transverse components, the inclusive cross 
section of single hard scattering reads 

^ps= 2 r^''^^i'2i)/(bi)'^t(^i.^'i) (7) 
xDI,{x\; e?)/(bi - b)dxidx[d^bid^b 



= D'h(xi\Ql)djk(xux\)Dl,(x\-Ql)dxidx\. 



These simplifying assumptions, though rather cus- 
tomary in the literature and quite convenient from a 
computational point of view, are not sufficiently justi- 
fied and are under the revision now lfl9[|20[|2lll27ll2il. 
However, the starting cross section formula dl] was 
found (derived) in many works (see, e.g., Refs. lliol 
U. 19. S 21]) in the momentum representation using 
the light-cone variables and the same approximations as 
thouse applied to the processes with a single hard scat- 
tering. 



Nevertheless, we restrict ourselves to this simple 
form (|5]l regarding it as the first estimation of the contri- 
bution from the double parton scattering to the inclusive 
double heavy meson production. The presence of the 
correlation term in the two-parton distributions results 
in the decrease ll26ll27ll28ll of the effective cross section 
CTeff with the growth of the resolution scales Qi and Q2, 
while the dependence of osff on the total energy at fixed 
scales is rather weak ll27ll . Thus, in fact, we obtain the 
minimal estimate of the contribution under considera- 
tion. 

The CDF and DO measurements give (Teff - 15 mb, 
which is roughly 20% of the total (elastic + inelastic) 
pp cross section at the Tevatron energy. We will use 
this value in our further estimations. 

Let us start from the double J/ 1// production, since the 
LHCb Collaboration has recently reported a first mea- 
surement |33] of this process 



= 5.6+1.1 ± 1.2 nb 



(8) 



with both y/i/r's in the rapidity region 2 < y < 4.5 
and with the transverse momentum pj"^ < 10 GeV/c 
in proton-proton collsions at a center-of-mass energy of 
yfs -7 Tev. Earlier this Collaboration has already mea- 
sured II34II the single inclusive J/tfr production cross sec- 
tion within the same kinematical cuts as above 



cr^ot = 7.65 ± 0.19 ± 1.10|'^f ^h. 



SPS 



(9) 



Using Eq. (|5]l we obtain immediately a simple esti- 
mation of the contribution from the double parton scat- 
tering at the same kinematical conditions 



DPS 



^ " SPS SPS 
2 CTeff 



^ 2.0 nb. 



(10) 



This value is quite comparable with the cross section 
through the "standard" mechanism of the double Jj^ 
yield |l3l|l 



cr^(f/^ = 4.15 nb, 



(11) 



and the theoretical prediction for the contribution from 
both scattering modes to the cross section 



" SPS " DPS 



6.15 nb. 



(12) 



is very close to the experimentally observed cross sec- 
tion (jSJ of double y/i/r production. 

A method to measure the double parton scattering at 
the LHCb using leptonic final states from the decay of 
two promt J 141 mesons is disccused in Ref. fy^. It 
is worth mentioning on the other hand that the predic- 
tions on the double // production are very sensitive to 



2 



the choice of the renormalization scale (because of the 
(9(aj) dependence of the cr^J^y^^'^ cross section), and so, 
the LHCb experimental results can also be accommo- 
dated by the SPS mechanism alone (see below). 

An even better evidence for the double parton scatter- 
ing process can be found in the production of Xc pairs. 
The production of P-wave states is suppressed relative 
to the production of 5 -wave states because of the hier- 
archy of the wave functions \'Rji^(Q)\^ » l*??^ (0)p/m^ 

leading to the inequality cr'^Jp^^^'^ » ■ Indeed, as 
one can learn from Fig. 8 in Ref. [3(ji], the inclusive 
double Xc production is suppressed in comparison with 
the inclusive double J/ production by more than two 
orders of magnitude. 

At the same time, the inclusive production of single 
J /iff andxc states shows nearly the same rates. The lat- 
ter property is supported by both theoretical llssl and 
experimental 11361 1371 13811 results. The reason comes 
from the fact that the Xc mesons are produced in a di- 
rect 2 — > 1 gluon-gluon fusion g + g ^ Xc, while 
the J/ if) mesons are produced in a 2 — > 2 subprocess 
g + g ^ J/ifr + g, where an additional final state gluon is 
required by the color and charge parity conservation. As 
a consequence, the invariant mass of the produced sys- 
tem is typically much higher in the J case than in the Xc 
case. (Besides that, the structure of the matrix element 
is such that it vanishes when the coproduced gluon be- 
comes soft. This further suppresses the production of 
low-mass states.) 

Taken together, the suppression factors coming from 
the lower wave function on the Xc side, and from the 
higher final state mass and extra as coupling on the 
y/i/r side nearly compensate each other making the in- 
clusive production cross sections comparable in size: 
cr^pj ^ cTsps- ^ consequence, we get cr^'ps - cr^'^''^'^ 



and cr^'pj » ai'^' 



DPS 



gP5 . Thus, if observed, the production of 
a XcXc pair would yield a clear and unambiguous indi- 
cation of the double parton scattering process. The need 
in detecting the decay photon ;(fc — » J/i// + j brings ser- 
tain difficulties in the experimental procedure, but the 
task seems still feasible as the production cross section 
is not small. 

Another tempting possibility is to consider the simul- 
taneous production of J/i// andxc- In the SPS mode, this 
process is forbidden at the leading order (LO) by the 
charge parity conservation but is possible at the next- 
to-leading order (NLO), g + g — > J/i(/ + Xc + g- The 
corresponding cross section is then suppressed by one 
extra power of as and by the Xc wave function. Al- 
ternatively, it can proceed via the soft final-state gluon 
radiation (the so called color octet model). The estima- 



tions of the cross section are then model dependent and 
rather uncertain, but even with the largest acceptable 
values for the color octet matrix elements one arrives at 
a suppression factor of about two orders of magnitude 
llsoll . For the DPS mode we still expect no suppression, 

J/^X. ^ J/'/'J/'l' 
DPS " DPS ■ 

Yet another interesting process is the production of 
particles from different flavor families, say, J/i// and 
Y mesons. Once again, this process is not possible at 
the leading order in the SPS mode and only can oc- 
cur either at the NNLO (next-to-next-to-leading order) 
0{a^s)' oi" via the color-octet transitions, or by means 
of the production and decay of P-wave mesons (i.e., 
g + g Xc + Xb followed by Xc Jl^ + 7 and 
Xb ^ + y)- So, the SPS mode is always suppressed: 
either by the extra powers of a,, or by the color-octet 
matrix elements, or by the P-state wave function, and 
the DPS mode becomes the absolutely dominant one: 

" DPS ^ " SPS ■ 

Now, to be more precise, we will derive some nu- 
merical predictions. In doing so, we rely upon per- 
turbative OCD and nonrelativistic bound state formal- 
ism 1 3^ 431 with only the color-singlet channels taken 



into consideration. Also, we accept the A;,-factorization 
ansatz ll4ll l42i l43ll for the parton model. The compu- 
tational technique is explained in every detail in Ref. 
113511, and the parameter setting is as follows. The me- 
son masses are taken from the Particle Data Book ll44ll . 
and the heavy quark masses are set equal to one half of 
the respective meson masses; the radial wave functions 
of J I if/ and T mesons are supposed to be known from 
their leptonic decay widths 144] and are set to |'??y/^(0)p 
= 0.8 GeV^ and |'??t(0)P = 6.48 GeV^; the wave func- 
tions of the P-states are taken from the potential model 
I??; (0)p = 0.075 GeV^ and \%^m^ = 1.44 GeV^; 
the renormalization scale in the strong coupling asijJ^) 
is chosen as the meson transverse mass /i^ = + pj; 
and we use the AO parameterization from Ref. 1461 for 
the unintegrated gluon density. 

In the present note we will restrict ourselves to the 
conditions of the LHCb experiment, since the Collabo- 
ration has already recorded the production of J/i// pairs. 
Predictions for other experimental conditions can be 
made in an essentially similar way. 

Within the theoretical model described above, we get 
for the direct inclusive y/f/r production 



(Tgpg (direct) = 7.1 yub, 
and for the Xc mesons 



£7^P5 = 1.5 fih, cr^pj = 5.1 fih. 



(13) 



(14) 



3 



After multiplying these numbers by appropriate 
branching ratios [44] Brixd^J/ifr+j) = 35% and 
Br(xc2-^J l^+y) - 20% and summing the direct and in- 
direct contributions together we get for the prompt J I tp 
yield 



Jl>l' 

SPS 



CTg^,'^ (direct) + a-i!^^{fxomxc) 



= 7.1/ib + 1.6;ub = 8.7//b. 



(15) 



This result is in reasonable agreement with the exper- 
imental measurement (9), thus giving support to our 
theoretical model. Quite similarly, we get for the bb 
mesons 



(Tsps (direct) = 140 nb, 



and 



"SPS 



18 nb, 



"SPS 



91 nb. 



Then one can easily obtain for the DPS mode 
cr^''p^/^^= 1.7nb, 
crQ'p^/^^(both from Xc) = 0.9 nb, 
o-q'p^^ = 0.07 nb. 



(16) 

(17) 

(18) 
(19) 
(20) 



The reader can continue deriving predictions for other 
DPS combinations. 

To calculate the background contribution cr^Jj!^^'^ we 
use the code developed in |47] and extended now i48ll 
to the A;,-factorization approach: 



u gpj t nu. 



(21) 



Variations in the renormalization scale yu^ within a fac- 
tor of 2 around the default value yu^ - i/4 make a fac- 
tor of 1.6 increasing or decreasing effect on the total 
production rate. Employing some different parametriza- 
tions for the unintegrated gluon densities (A+ or A- sets 
from Ref. (A6\) also changes the predicted cross section 
by a factor of 1.6 up or down. Our central prediction 
(21) is in reasonable agreement with the data (8). 

The proportion between the visible SPS and DPS 
contributions can, in principle, depend on the experi- 
mental cuts on the Jjili transverse momentum. However, 
in the particular case which we are considering here, the 
LHCb Collaboration refers to no cuts on priJ/i//)- In 
fact, there are some soft restrictions on the momenta of 
the decay muons, pril-i) > 600 MeV, but they are taken 
into account as corrections to the acceptance. The final 
results reported by the Collaboration to compare with 
are the acceptance-connected ones. 

It is also worth noting that even in the general case the 
sensitivity of the ratio o-ups/osps to the pr cuts is rather 



weak, because the DPS and SPS contributions show the 
same pr behavior This is explained in detail in Ref. 
|48]. Irrespective of the particular properties of the sub- 
process matrix element, the pr of the final state is domi- 
nated by the transverse momentum of the initial gluons, 
and the individual J/ 1// spectra behave as 1 / pj in both 
SPS and DPS modes. Moreover, the momenta of the 
two J/i// mesons are not correlated. The latter is evi- 
dent in the DPS case and was not a priori evident in the 
SPS case, but turned out to be true (Fig.5 in Ref. iH]). 
So, the SPS and DPS event topologies are rather similar 
to each other and can hardly be distinguished from one 
another 

Our calculations agree with the observation made in 
Ref. 1321 that the effects of initial parton radiation 
(that are automatically present in the A;rfactorization ap- 
proach) destroy the original back-to-back kine- 
matics completely washing out the azimuthal correla- 
tions. One can potentially distinguish the SPS and DPS 
modes with rapidity correlations, but we anyway find 
that looking at some other meson species is more in- 
dicative. In particular, the production of XcXc, JI^Xc 
or Jjijjy pairs is totally dominated by the DPS mecha- 
nism because the SPS mechanism is suppressed for the 
reasons given earlier 

Summing up, we conclude that the processes with 
pairs of heavy quarkonia in the final state {J/iJ/J/il/, 
XcXc, JI4'Xc^ J/^'^ ) can serve as precise probes of the 
double parton scattering at the LHC and can stimulate 
important steps towards understanding the multiparticle 
QCD dynamics. 
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